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Abstract 

Constraints on the general 2HDM ("Model II") are obtained from 
the existing {g — 2)« data including limits on Higgs bosons masses 
from LEP I data. We consider separately two cases: with a light 
scalar h and with a light pseudoscalar A, assuming Mh + Ma > Mz- 
The charged Higgs contribution is also included. It is found that 
already the present (g — 2)^ data improve limits obtained recently 
by ALEPH collaboration on tan (3 for the mass of the pseudoscalar 
below ~ 2 GeV. The improvement in the accuracy by factor 20 in 
the forthcoming E821 experiment may lead to more stringent, than 
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provided by ALEPH group, limits up to Ma ~ 30 GeV if the mass 
difference between h and A is ~ Mz- Similar results should hold for 
a light scalar scenario as well. 

1 Status of 2HDM. 
1.1 Introduction. 

The mechanism of spontaneous symmetry breaking proposed as the source 
of mass for the gauge and fermion fields in the Standard Model (SM) leads 
to a neutral scalar particle, the minimal Higgs boson. According to the LEP 
I data, based on the Bjorken process e + e~ — > HZ*, it should be heavier than 
66 GeV0, also the MSSM neutral Higgs particles have been constrained by 
LEP1 data to be heavier than ~ 45 GeV 0, |3], [p. The general two Higgs 
doublet model (2HDM) may yet accommodate a very light ( ^ 45 GeV) 
neutral scalar h or a pseudoscalar A as long as M h + Ma ^ M z @ . 

The minimal extension of the Standard Model is to include a second Higgs 
doublet to the symmetry breaking mechanism. In two Higgs doublet models 
the observed Higgs sector is enlarged to five scalars: two neutral Higgs scalars 
(with masses Mh and Mh for heavier and lighter particle, respectively), one 
neutral pseudoscalar (Ma), and a pair of charged Higgses (M H + and Mh-). 
The neutral Higgs scalar couplings to quarks, charged leptons and gauge 
bosons are modified with respect to analogous couplings in SM by factors that 
depend on additional parameters : tan/3, which is the ratio of the vacuum 
expectation values of the Higgs doublets t^/fi, and the mixing angle in the 
neutral Higgs sector a. Further, new couplings appear, e.g. Zh(H)A and 
ZH+H-. 

In this paper we will focus on the appealing version of the models with 
two doublets (" Model II" ) where one Higgs doublet with vacuum expectation 
value t>2 couples only to the "up" components of fermion doublets while the 
other one couples to the "down" components ||. In particular, fermions cou- 
ple to the pseudoscalar A with a strength proportional to (tan/3) ±1 whereas 
the coupling of the fermions to the scalar h goes as ±(sina/ cos/?) ±1 , where 
the sign ± corresponds to isospin =Fl/2 components. In such model FCNC 
processes are absent and the p parameter retains its SM value at the tree 
level. Note that in such scenario the large ratio v^/vx ~ mtop/rrib ^ 1 is 
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naturally expected. 

The well known supersymmetric model (MSSM) belongs to this class. In 
MSSM the relations among the parameters required by the supersymmetry 
appear, leaving only two parameters free (at the tree level) e.g. M4 and tan (5. 
In general case, which we call the general 2 Higgs Doublet Model (2HDM), 
masses and parameters a and (3 are not constrained by the model. Therefore 
the same experimental data may lead to very distinct consequences depending 
on which version of two Higgs doublet extension of SM, supersymmetric or 
nonsupersymmetric, is considered. 

1.2 Present constraints on 2HDM from LEP I. 

Important constraints on the parameters of two Higgs doublet extensions of 
SM were obtained in the precision measurements at LEP I. The current mass 
limit on charged Higgs boson M H ±= 44 GeV/c was obtained at LEP I [|J 
from process Z — > H + H~, which is independent on the parameters a and 
(3. (Note that in the MSSM version one expect M H ± > Mjy). For neutral 
Higgs particles h and A there are two main and complementary sources 
of information at LEP I. One is the Bjorken processes Z Z*h which 
constrains g\zz ~ sin 2 (a ~ (3)- The second process is Z — > hA, constraining 
the g ZhA ~ cos 2 (a — 0) for M h + Ma~M z []. This Higgs pair production 
contribution depends also on the masses M h , Ma and Mz- 

Results on sin 2 (a — (3) and cos 2 (a — (3) can be translated into the limits 
on neutral Higgs bosons masses Mh and Ma- In the MSSM, due to relations 
among parameters, the above data allow to draw limits for the masses of 
individual particles: Mh > 45 GeV for any tan (3 and Ma > 45 GeV for 
tan/? >1 @, In the general 2HDM the implications are quite different, 
here the large portion of the (M^M^) plane, where both masses are in the 
range between and ~50 GeV, is excluded M. 

The third basic at LEP I process in search of a neutral Higgs particle 
is the Yukawa process, i.e. the bremsstrahlung production of the neutral 
Higgs boson h(A) from the heavy fermion, e + e~ —> ffh(A), where / means 
here b quark or r lepton. This process plays a very important role since 
it constrains the production of a very light pseudoscalar even if the pair 

1 The off shell production could also be included, e.g. as in j|. 
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production is forbidden kinematically |, i.e. for Mh + Ma > Mz- It allows 
also to look for a light scalar, being an additional and in case of a = /3 the 
most important, source of information. The importance of this process was 
stressed in many papers || |7j, the recent discussion of the potential of the 
Yukawa process is presented in Ref.|8|]. 

New analysis of the Yukawa process by ALEPH collaboration || led to 
the exclusion plot (95%) on the tan/3 versus the pseudoscalar mass, Ma- 
(Analysis by L3 collaboration is also in progress fL0| .) It happened that 



obtained limits are rather weak 0, allowing for the existence of a light A with 
mass below 10 GeV with tan/3 = 20-30 , for M A =40 GeV tan/3 till 100 is 
allowed ! For mass above 10 GeV, similar exclusion limits should in principle 
hold also for a scalar h with the replacement tan/3 — > sin a/ cos (3. Larger 
differences one would expect however in region of lower mass, where the 
production rate for the scalar is considerably larger than for the pseudoscalar 
and therefore more stringent limits should be obtained ||. 

1.3 The 2HDM with a light Higgs particle. 

In light of the above results from precision experiments at LEP I there is 
still a possibility of the existence of one light neutral Higgs particle with 
mass below ~ 40-50 GeV. As far as other experimental data are concerned, 
especially these from low energy measurements, they do not contradict this 
possibility as they cover only part of the parameter space of 2HDM, moreover 
some of them like the Wilczek process have large theoretical uncertainties 
both due to the QCD and relativistic corrections |L1, || (see also discussion 
in 



12, 13, 14 



In following we will study the 2HDM assuming that one light Higgs parti- 
cle may exist. Moreover we will assume according to LEP I data the following 
mass relation between the lightest neutral Higgs particles: M^ + Ma > Mz- 
We specify the model further by choosing particular values for the param- 
eters a and (3 within the present limits from LEP I. Since sin(a — ft) 2 was 
found [0, H to be smaller than 0.1 for the 0~M/i~ 50 GeV, and even below 
0.01 for a lighter scalar, we simply take a = (3. This assumption leads to 
equal in strengths of the coupling of fermions to scalars and pseudoscalars. 



2 neglecting the off shell production 

3 Note, that the obtained limits are much weaker than the limits estimated in Rcf. 
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For the scenario with large tan/3 ~ Oimtjm^) large enhancement in the cou- 
pling of both h and A bosons to the down-type quarks and charged leptons 
is expected. 

As we described above the existing limits from LEP I for a light neutral 
Higgs scalar /pseudoscalar boson in 2HDM are rather weak. Therefore it is 
extremely important to check if more stringent limits can be obtained from 
other measurements. The possible effect due to such a light Higgs particle 
at ep collider HERA and at low energy LC, as well in heavy ion collisions at 



HERA and LHC are discussed elsewhere [12|, [L5|, [16], [L7| • 

In this paper we study in details the limits on the 2HDM from the pre- 
cision {g — 2) experiment for muon[]. We have studied this problem earlier 
using the simple approach constraining the individual contributions h (or A) 
(see Ref.[|T3l, [14], [jjj and also 0). In the present analysis we take into account 
the full contribution from 2HDM, i.e. exchanges of h, A and bosons in- 
corporating the present constraints on Higgs bosons masses from LEP I. We 
study here the present data on (g — 2) M flTj| as well as the potential of the 
future E821 experiment JT9[ with the accuracy expected to be more than 20 
times better than the present one. 



2 Constraints on the parameters of 2HDM 
from (g — 2) for the muon. 

2.1 The room for new physics. 

The present experimental data limits on (g — 2) for muon averaged over the 
sign of the muon electric charge is given by ||21[| : 



a exp _ (9 ^ = 1 lg5 g23 (g^ . 1Q -9 

The quantity within parenthesis, cr exp , refers to the uncertainty in the last 
digit. The expected high-precision E821 Brookhaven experiment has design 
sensitivity of c^™ = 4 ■ 10~ 10 (later even 1-2 -10~ 10 |33]]) instead of the above 



84 • 10 10 . It is of great importance to reach this accuracy in the theoretical 
analysis. 

4 Measurement of (g — 2) of the electron does not give useful restriction due to small 
mass of the electron (see also p5[). 
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The theoretical prediction of the Standard Model for this quantity consist 
of the QED, hadronic and EW contribution: 



,SAf 



OED + Kad + a EW_ 

H 1 p 



The recent SM calculations of a 



M are based on the QED results from |22|, g3 

and 



25, 24, 26 



33 



31 and 



hadronic contribution obtained in 
the EW results from |33|, |32j. The uncertainties of these contributions differ 
among themselves considerably (see e.g. discussion below and in Ref. [20 
13 



|27j|). The main discrepancy is observed for the hadronic contribution, 
therefore we will consider here two cases : case A f\ with relatively small 
error in the hadronic part and the case B Q with the 2 times larger error in 
the hadronic part. (We adopt here the notation from pOf but one should 
be aware that the numbers used in this analysis differ slightly from Ref. [^(J . 
Basically our case A corresponds to case A from f20fl , whereas case B based 
on the fl33fl.) 



case A [in KT 9 ] B [in 1(T 9 ] 

QED 1 165 847.06 (0.02) 1 165 847.06 (0.02) 

had 69.70 (0.76) 68.82 (1.54) 

EW 1.51 (0.04) 1.51 (0.04) 

Tot 1 165 918.27 (0.76) 1 165 917.39 (1.54) 



Note that the hadronic contribution error dominates the total error for 
the SM predictions and will influence strongly the comparison with the new 
precision data from E821, what will be discussed later. 

The room for a new physics is given basically by the difference between 
the experimental data and theoretical SM prediction: a e * p — a^ M = 5a^. In 
the following we will assume that it is consistent to use this difference 5a^ as 
the indication for the contribution due to Higgs particle (s) of a beyond SM 
origin, like 2HDM, although in the calculation of a^ w the (SM) Higgs scalar 
contribution is included. This assumption based on the observation that two 
EW calculations: the EW result from Ref. |33j based on M#j 99S =250 GeV and 



5 Refs. 
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6 Refs. 


22 





24 
31 



30, B3 
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the one from Ref.|2| with M Higgs > 5 GeV, differ only by ~ 0.02(0.06) • 10~ 9 . 



Below the difference Sa^a) for these two cases: A and B, is presented 
together with the error a, obtained by adding in quadrature the experimental 
and theoretical errors: i.e. a = J (a. 



exp 



°ix P + o- 



2 1 

had I 



case A [in 10~ 9 ] B [in 10~ 9 ] 

Sa^ja) 4.73(8.43) 5.61(8.54) 

lim(95%) -11.79 < Sa^ < 21.25 - 11.13 < <5a M < 22.35 

lim±(95%) -13.46 < 5a^ < 19.94 - 13.71 < <?a M < 20.84 



One can see that at one sigma level the difference Sa^ can be positive or 
negative. For that beyond SM scenarios in which both positive and nega- 
tive may appear, the 95% confidence level (C.L.) bound can be calcu- 
lated straightforward (above denoted by lim(95%)). For the model where 
the contribution of only one sign is physically accessible, the other sign be- 
ing unphysical, (i.e. positive or negative 5a^) the 95%C.L. limits should be 
calculated in different way ||21||. These limits calculated separately for the 



positive and for the negative contributions are denoted above by Um±(95%). 

We found that this latter approach leads to the sizeable shift in the lower 
(/_) and upper limits with respect to the standard (95%) limits (by -1.3 
•10~ 9 up to -2.6 -10~ 9 ). That means that the possible negative contribution 
becomes larger whereas the positive becomes smaller when lim± method is 
used. The differences between theoretical predictions (case A and B) for fixed 
method of calculating confidence level seem to be smaller that this effect. All 
these effects may be important in future analysis. 



2.2 2HDM contribution to (g - 2) M . 

As we mentioned above the difference between experimental and theoretical 
value for the anomalous magnetic moment for muon we ascribe to the 2HDM 
contribution, so in order to constrain the parameter space of the model we 
take5a^ = a$ HDM \ 

7 The contribution due to the Higgs scalar for the Mniggs—5 GeV was found to be of 
the level of 1 -lCT 11 Q. 
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To aff HDM ' a scalar h (ajj), pseudoscalar A (a^) and the charged Higgs 
boson if ± (ajj) contribute. The relevant formulae can be found in the Ap- 
pendix. Each term (A = h, A or H^) disappears in the limit of large 
mass, at small mass the contribution reaches its maximum (or minimum if 
negative) value. The scalar contribution a^(M fe ) is positive whereas the pseu- 
doscalar boson ci^(Ma) gives negative contribution, also the charged Higgs 
boson contribution is negative. Note that since the mass of H ± is above 44 
GeV (LEP I limit), its small contribution can show up only if the sum of h 
and A contributions is small. 

We calculate the a@ HDM > minimalizing the cross section in order to put 
limit on the maximum allowed tan/5. Therefore we take M°=44 GeV and 
for mass of neutral Higgs bosons we assume the relation Mh + Ma > M z as 
discussed before. 

We study two scenarios: 

• a) pseudoscalar A is light, the scalar mass M h > M z — Ma and we 
calculate the total 2HDM contribution as follows: 

af DM) (M A ) = o*{M A ) + aJ(M° = M z - M A ) + aJ(M°) (la) 

• b) scalar h is light, the pseudoscalar mass is Ma > M z — M h and: 

a (2HDM ){Mh) = a H (Mh) + a A (M A = Mz _ + a ± (M ± y ^ 

Due to opposite signs there appear a cancellation in both scenarios be- 
tween scalar and pseudoscalar contributions, especially strong for Ma ~ 
in (la) or M h ~ M° A in (lb). Note that the total 2HDM contribution is for 
the scenario ajnegative, whereas for the scenario b)- positive. Therefore we 
have to include this fact calculating the 95% C.L. bounds of af^ HDM \ This 
leads us to the limits Zim±(95%) described previously in Sec. 2.1. 

Present data. 

8 the case with mass equal 600 GeV was also studied, but in the range of masses of 
neutral Higgs bosons below 40 GeV discussed here, the influence of so heavy charged 
Higgs boson is negligible, see below. 
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Since the case A (Sec. 2.1) gives more stringent lim±(95%) constraints, 
for both positive and negative contributions, only this case will be used in 
the further analysis. So, we have for the considered 2HDM scenarios allowed 
following ranges: 

a) -13.46-HT 9 < a% BDM ) < (2a) 

b) < af HDM ^ < 19.94-10- 9 . (26) 

Using Eqs.2a(2b) the 95% C.L. exclusion plots for tan/5 versus Ma (Mh) 
is calculated for a light pseudoscalar (light scalar) scenario. The results are 
presented in Fig.l (2). The regions above curves are excluded. 

The results based on the current data for (g — 2) M are presented by upper 
curves. The total contribution of the 2HDM (solid line) as well the simple 
approach where only pseudoscalar (scalar) contribution (short-dashed line) 
is included are presented in Fig.l and 2. They lead to similar limits up to 
mass ~ 20 GeV. Note that the simple approach for pseudoscalar (scalar) 
corresponds to the limits of negligible contribution of scalar (pseudoscalar), 
i.e. to the large difference in mass between h and A (above 150 GeV or so). 

The charged Higgs boson contributes very little (in Figs.l and 2 the dif- 
ference between a solid line and the long-dashed line, calculated without the 
charged Higgs boson term) being visible only for masses of h and A above 40 
GeV, where the cancellation between scalar and pseudoscalar becomes very 
strong (see small figures in Figs.l, 2 for details). 

Interestingly, the present (g — 2)^ data can accommodate large value of 
tan/3 (20 or more) for the Higgs boson masses equal or larger than 2 GeV 



(see also the discussion in Ref. [0, [13], (To)). 

Forthcoming data. 

Since the dominate uncertainty in 5a^ (Sec. 2.1) is due to the experimental 
error, the role of the forthcoming E821 experiment is crucial in testing the SM 
or probing a new physics. We discuss now the potential of this measurement 
for the constraining the 2HDM. 

The future accuracy of the (g — 2) M experiment is expected to be a™p ~ 
0.4 ■ 10~ 9 or better, and one may in principle calculate the expected error for 
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the difference 8a n ™ . Assuming no progress in theoretical SM calculation the 
above experimental accuracy will lead to a new = 0.86 and 1.6, for case A and 
B, respectively. 

However, especially for the case A discussed in Sec. 2.1, one could in prin- 
ciple expect such an improvement in the calculation of the hadronic con- 
tribution that the total uncertainty in Sa^ will be basically due to the 
experimental error. We will take this point of view and in particular we will 
assume that the accessible range for a$ HDM \ given presently by Eq.2a(2b) 
for a light pseudoscalar (scalar), would be smaller by factor 20. In this ap- 
proach the central value for the difference da^ is shifted by the same factor 
20 to lower value 0.24 -10~ 9 (from 4.73 -10~ 9 ). So, we consider the following 
option (in 10~ 9 ): 

5a n ™ = 0.24, and lii< cw (95%) : -0.69 < 5a^ < 1.00 

or separately for two scenarios: 

a') -0.69 -10- 9 <a^ HDM ^ <0 (3a) 

b') < a^ HDM) < 1.00- 10" 9 (36) 

As expected much stringent limits on tan/3 are obtained in this case for 
pseudoscalar (scalar) - see Fig.l (2) lower curves. 

2.3 Discussion and comparison with ALEPH limits. 

The comparison of the obtained 95%C.L. exclusion plots based on the current 
and future accuracy for the (g — 2) M data with the latest results from ALEPH 
collaboration based on the Yukawa process |J] is presented in Fig. 3. In this 
figure a simple approach prediction for a light h (solid line) or a light A (long- 
dashed line) are shown for current and future (g — 2)^ measurement, upper 
and lower curves, respectively. 

9 The improvement in the ongoing experiments at low energy in expected as well. 
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The results from ALEPH collaboration are presented by dotted line, only 
for the pseudoscalar case. At very low mass, below say 2 GeV, the present 
(g — 2)^ data leads to stronger limits on tan f3 than the Yukawa process. 

Similar results from Yukawa process as presented by dotted line should 
also hold for the case of a light scalar for M h above 10 GeV. For lower mass 
both Yukawa and (g — 2)^ data may lead to comparable exclusions plots. 

Note that the predictions based on the simple approach can be justified 
if the mass difference between h and A is bigger than, say 150 GeV. In this 
case the improved (g — 2) M data for muon will lead to more stringent limits 
over the whole considered range of masses. If the mass difference between h 
and A is ~ M z the improvement may be obtained up to ~ 30 GeV. 

3 Conclusion 

We studied the room for new physics in the current muon anomalous mag- 
netic moment data based on the newest theoretical calculations of the SM 
prediction a^ M \ 

The difference between experimental and theoretical value for the anoma- 
lous magnetic moment for muon was ascribed to the 2HDM contribution, so 
we took 5a^ = a^ HDM \ From evaluated 95% C.L. bounds for the positive 
or negative value of 

a (2HDM) we 

derived constraints on the general 2HDM 

("Model II"). 

We studied the total 2HDM contribution aff HDM ^ separately in two cases: 
with a light scalar h (M h < 45) and with a light pseudoscalar A (M h < 45), 
assuming according to the LEP I data that M h + Ma > M z . A light scalar 
scenario leads to the positive, whereas the one with a light pseudoscalar to 
the negative a^ HDM \ Large cancellation occure when both mass of neutral 
Higgs bosons are nearly degenerated M h ~ M A ~ M z /2. In this case the 
charged Higgs boson contribution may show up in the range of M h (A) above 
40 GeV, if we take the lower LEP I limit for M ± =44 GeV. 

We found that the contribution due to a light A becomes larger whereas 
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this one due to a light h becomes smaller by ~ 1.5 or even 2.6 -10 -9 if the 
proper evaluation of 95% CL. is introduced. The differences between theo- 
retical predictions (mainly due to the uncertainty in the hadronic corrections) 
on the level of 0.25 to 0.9 -10~ 9 are smaller that this effect. Note also that 
the the theoretical uncertainty influences less the lower limit, relevant for 
the pseudoscalar contribution, than the upper one, where a scalar contribu- 
tion is expected. All these effects may be important for the future accuracy 
measurement of (g — 2) for muon. 

It was found that already the present data improve limits obtained 

recently by ALEPH collaboration on tan f3 for low mass of the pseudoscalar: 
Ma < 2 GeV. The future improvement in the accuracy by factor 20 in the 
forthcoming E821 experiment may lead to more stringent limits than pro- 
vided by ALEPH group up to M^= 30 GeV or higher for a larger mass 
difference between scalar and pseudoscalar. Similar results should hold for a 
2HDM with a light scalar. 
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5 Appendix 



Theoretical predictions of an anomalous magnetic moment of muon from 
triangle vertex with the exchange of the Higgs boson A (= h, A or H^) are 



given in Ref . |36|, |35 



/a 



57T 
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where the coupling constant /a is given by 

I A = — tan 3, 

2 M w 11 

(in case of scalar tan/5 should be replace by the sin(a)/ cos(/3), but in our 
model we put a = (3, so the coupling constant is universal ). 

The integral Lu(A) is for the neutral Higgs contribution given by: 



/' 

Jo 



" 1 " <,X x* + (l-x)(M h/A /m,) 2 ' 
with 



Q h (x) = x 2 (2-x) (1) 
Q A {x) = -x\ (2) 



The charged Higgs particle exchange is described by: 

L ± = j l dx - X{1 ~ X 1 . 
Jo x + (M ± /m^) 2 - 1 
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Exclusion plot, 95% C.L: light pseudoscalar 
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Figure 1: The 95% exclusion plot, based on lim±, from the (g — 2)^ data 
for pseudoscalar. The present limits (Eq.2a) (upper curves) and for the im- 
proved measurement (Eq.Sa) (lower curves) are presented. The total 2HDM 
(Eq.la) contribution is represented by the solid lines, the short-dashed lines 
correspond to the simple approach where only pseudoscalar contributes. The 
long-dashed lines correspond to the case where the charged Higgs particle con- 
tribution is neglected; for details see small figure, where the mass range above 
44 GeV is displayed. 
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Exclusion plot, 95% C.L.: light scalar 
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Figure 2: The 95% exclusion plot, based on lim±, from the (g — 2)^ data 
for scalar. The present limits (Eq.2b) (upper curves) and the improved mea- 
surement (Eq.Sb) (lower curves) are presented. The total 2HDM (Eq.lb) 
contribution is represented by the solid lines, the short-dashed lines corre- 
spond to the simple approach where only scalar contributes. The long-dashed 
lines correspond to the case where the charged Higgs particle contribution is 
neglected; for details see small figure, where the mass range above 44 GeV is 
displayed. 
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Figure 3: The 95% exclusion plot, based on lim±, from the present (upper 
curves) and future (lower curves) (g — 2)^ data. The simple approach results 
for a light h (solid line), and a light A (long-dashed line) with the new ALEPH 
analysis for pseudoscalar (dotted line) are presented. 



19 



